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Abstract 

q is the average squared transverse momentum transfer per unit length to a high-energy particle 
traversing a QCD medium such as a quark-gluon plasma. We find the (UV-regulated) value of q to 
leading order in the weak coupling limit, a s (T) <C 1- We then use this value to generalize previous 
analytic results on the gluon bremsstrahlung and pair production rates for massless high-energy 
particles in a weakly-coupled quark-gluon plasma, at next-to-leading logarithmic order. 



I. INTRODUCTION 



In a weakly-coupled quark-gluon plasma, energy loss of high-energy jets is dominated 
by gluon bremsstrahlung and pair production. One of the elements that goes into the 
calculation of gluon bremsstrahlung rates in the high-energy limit is the parameter q — the 
average squared transverse momentum transfer per unit length to a high-energy particle 
due to elastic collisions with plasma particles. In this paper, we find the (UV-regulated) 
value of q to leading order in the weak coupling limit. We then use it to generalize previous 
analytic results on the gluon bremsstrahlung and pair production rates at leading order in 
coupling and next-to-leading order in inverse powers of \n(E/T), where E is the energy of a 
jet particle. 



A. The parameter q 

We will define the parameter q by 

qW= f <?q ± ^ql, (1.1) 

Jq ± <A a Q± 

where T e \ is the rate for elastic collisions with plasma particles, q± is the transverse momen- 
tum transfer in a single such collision, and A is an ultraviolet cut-off whose physical scale 
will be set later when we discuss gluon bremsstrahlung. The value of q is proportional to a 
factor of the quadratic Casimir Cr for the color representation R of the high-energy particle. 
It will be convenient to factor out this dependence by defining q by 

q = C R l (1.2) 

The differential elastic cross-section has the known limiting forms 1 

qKql + ml)' ± 



dT e \ _ Cr 

X ^ 



d 2 q± (2n)< 



where m D is the Debye mass and Af is the density of plasma particles times color group 
factors (details below). 

If we wanted q(A) for A <C T, we could simply use the first case of (1.3) to obtain 2 

m -L^ffA- 2aT< in (4) (mD * A<<n {iA) 

However, as we shall discuss momentarily, bremsstrahlung for sufficiently high energy parti- 
cles depends on if (A) for A>T. For this, we will need to find dT e \/d 2 q± for q± ~ T as well. 



1 The simple formula for the small q± case is taken from Aurenche, Gelis, and Zakaret [2]. The dT e \/d 2 q± 
presented here corresponds to the notation g 2 Cn A(q±) / (2tt) 2 of Refs. [6, 11]. See Appendix A of Rcf. 
[13]. 

2 For (1.4), see also Eq. (13) of Ref. [14] and the relation to Rcf. [15] discussed after Eq. (61) of Ref. [14]. 
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Our final result for q(A) in this case is given in Sec. II below, and the derivation follows in 
Sec. III. 

Finding q(X) for A ^> T is an interesting problem of principle: It's always good to 
understand how to get precise answers to questions in the limit of weak coupling. But it's 
important to note that, even within the context of the weak coupling limit, getting it right 
will not have a large effect on the answer. As noted in Ref. [11], the numerators g 2 m^T 
and g A N in (1.3) differ by only about 15% for 3-flavor QCD. So one would not make a large 
error by using the simple formula (1.4) even when A>T. Nonetheless, in the present work 
we aim to find the exact weak-coupling answer for this limit. 

B. Gluon bremsstrahlung 

When very high-energy particles travel through a weakly-coupled quark-gluon plasma, 
the dominant energy loss mechanism is through hard bremsstrahlung or pair creation. 
Each time a high energy particle collides with a plasma particle, there is the potential 
for bremsstrahlung or pair creation. At high energies (parametrically E > T), the quan- 
tum mechanical duration (formation time) of a splitting process becomes larger than the 
mean free time between the underlying collisions with plasma particles, leading to coher- 
ence effects that reduce the bremsstrahlung or pair production rate. This is known as the 
Landau- Pomeranchuk-Migdal (LPM) effect. In this paper, we consider the calculation of 
such splitting rates for particles with large energies E ^> T, where T is the temperature of 
the plasma. We will consider the case of a time-independent, uniform, infinite quark-gluon 
plasma. In practice, that means that the plasma temperature should not vary significantly 
over the time and distance scales associated with the formation time of the splitting process. 

A calculation of these rates to leading-order in a s was performed by Jeon and Moore 
[3], based on the formalism of Arnold, Moore, and Yaffe (AMY) [4-6]. These calculations 
require numerical solutions of integral equations. It was noted earlier by Baier, Dokshitzer, 
Mueller, Peigne and Schiff (BDMPS) [7-10], however, that analytic solutions can be found 
if the logarithm of the energy is large, and they performed a leading-log analysis. Arnold 
and Dogan [11] recently extended this analysis to next-to- leading order in inverse powers of 
\n(E/T). The results of all of these analysis depend on the differential rate dT e \/d 2 q±_ for 
elastic scattering of the high energy particle off of a plasma particle, as a function of the 
transverse momentum exchange qj_. The full numerical calculations of rates to leading order 
in a s by Jeon and Moore [3] , and its approximation to next-to- leading log (NLL) by Arnold 
and Dogan [11], were carried out using the q± <C T case of (1.3). However, the relevant 
range of q± which contributes to the calculation grows as the energy E of the splitting parton 
grows. Roughly speaking, as the formation time gets longer, that time subsumes more and 
more elastic collisions, which means that the total momentum transfer Qj_ from the plasma 
to the splitting particle during the formation time grows larger. Q± grows parametrically 
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as 3 

Q± ~ (qE) 1/4 - (1-5) 
The total momentum transfer Q± is made up of many individual momentum transfers q±, 
which range in scale from m D to Q± itself. At high enough energy, the assumption q±_ < T 
therefore breaks down for the upper end of this range. Based on (1.5), the condition Q± <C T 
for assuming all q± <C T is 

£<— . (1.6) 
Q 

Using (1.4) and m D ~ gT, this condition is parametrically [11] 

E<< 7Wiiv < 17) 

In the current work, we will redo Arnold and Dogan's next-to-leading log calculation of hard 
bremsstrahlung and pair production, rewriting the answer in terms of q. We will then be 
able to use our result for q(A) with A 3> T to generalize the previous results to the new case 

E>> ¥WhY m 

Our result is given in Sec. II below, and the few modifications to the original derivation of 
Arnold and Dogan [11] are outlined in Appendix D. 



C. Weak coupling and mo <C T 



In this paper, we ruthlessly work in the weak coupling limit. In particular, the Debye mass 
mD is parametrically of order gT, and so we shall formally assume that m D <C T. In practice, 
however, one interest in weak coupling calculations is to see what they give if optimistically 
applied to realistic situations where the coupling is not terribly small. Unfortunately, the 
weak coupling result for for massless 3-flavor QCD gives, for example, tod ~ 2.4T 
at a s ~ 0.3. Exactly how terrible it is to treat m D <C T depends on the details of the 
calculation. In this paper, we simply explore the weak coupling limit and will not consider 
possible improvements that might be made by not assuming mo <C T. Recent progress 
on going beyond this approximation has been made by Caron-Huot [1], who pushes the 
calculation of dY e \/d 2 q^_ and q to next-to-leading order in coupling g. 



II. RESULTS 



A. Notation 



Throughout, we use the same notation for group factors as in Ref. [11]. For a color 
representation R, Cr is the quadratic Casimir, da is the dimension, and tn = d^Cn/dA is 



3 Eq. (1.5) is for the case where the gluon momentum fraction x ~ 1 and 1 — x ~ 1, which dominates 
energy loss. More generally, the scale of the formation time is determined by replacing E by the smallest 
of the final particle energy and bremsstrahlung gluon energy. The estimate (1.5) can be understood as 
Q\ ~ gtform, where t[ orm is the formation time, combined with the standard result that t[ orm ~ y/E/q 
for x ~ 1 — x ~ 1. See, for example, Sec. 3 of Ref. [8], and Ref. [12]. 
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the trace normalization. For QCD, 

Ca = 3, Cf = 3 , cIa = 8, g?f = 3, tA = 3, tp = \ (2.1) 

for gluons (A) and quarks (F). 

We will use S b and Sf to represent the number of spin and flavor degrees of freedom in 
plasma bosons and fermions, respectively, times the corresponding trace normalizations £r. 
In A f -flavor QCD, 

H b = it a = 6, S f = 4A f t F = 2JV f . (2.2) 

The Debye mass and the weighted density M appearing in (1.3) are given by the following 
formulas, which are cast in a form that will be useful later on: 

ml = [~ h ( + (2) + SfC_(2)] ^ = (1 + |iV f )/T 2 (2.3a) 

and 

N = [H b C + (3) + S f C_(3)] ^ = || (1 + |A f )T 3 . (2.3b) 
Here the functions (±{z) are the bosonic and fermionic versions of the Riemann ( function, 

^)-e4^' (2 - 4) 

fc=i 

and so 

C + («) = C(s), C-(a) = (l-2 1 -)C(a)- (2-5) 

Recall that C(2) = tt 2 /6. 



B. Result for q 

Our final result for q(A) for large cut-off A is 



with 



and 



|(A) = [S b 2 + (A) + S f J_(A)]: 



7T 



(A>T) 



X ± (A)^%2ih, 



A J+ = 



C±(2) ~ C±(3) 
2tt 



2tt 



In 



A 

m D 



+ AJ ± , 



— ) +|-7 E + ln2 



2tt' 



fc=l 



(2.6a) 

(2.6b) 
(2.6c) 

(2.6d) 



Above, 7 E is the Euler-Mascheroni constant. As we'll discuss later, the sums a± are related 
to a certain generalization of the ( function, but we are unaware of any way of writing them 
in terms of more usual mathematical constants. Their numerical values are 



a+ = 0.386043817389949- 



and 



cr_ = 0.011216764589789' 



(2.7) 
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The result just quoted ignored running of the coupling constant. If A is so large that g 2 {A.) 
is significantly different from g 2 {mY)) and g 2 (T), then one should make the replacements 



,4, ( A \ . 2/ A x 2/ , , ( A 



/In — ^(A)^(m D )ln — , (2.8) 
5 4 lnf— "I ^^(T)^ 2 (m D ) \n(—) , (2.9) 



\m D y \m D J 

following the discussion in Refs. [11, 13, 17], and replace the g 4 which do not multiply logs by 
g A {T) as these constants turn out to be determined by the q±~T range of the q± integration 
in (1.1). This compact prescription accounts for 1-loop running of the coupling constant, 
provided there are no vacuum mass thresholds between and A. The difference between 
g 2 {T) and g 2 (m D ) is of order g 4 ln(T/m D ) ~ g A ln(g _1 ) and so not significant in the weak 
coupling limit, but no harm is done in accounting for this particular higher-order correction 
to q. 

The replacement (2.8) has a finite limit as A — > oo: 

94] /M ^) (A = oo)> (2 . 10) 



\m D J -2{3 



where 

(11 C A - 4N { t F ) (33 - 2 N { 
48tt 2 " 48^ 



po — 2 — j^- 2 — y 2 - 11 ) 



is the one-loop coefficient of the /3 function for g 2 (fi) = [— (3q ln(/i 2 /A 2 )]" 1 . So, if the effect 
of running is included, q{oo) is finite and given by the above formulas. For applications to 
bremsstrahlung, however, we will generally want <f(A) for finite A. 



C. Result for gluon bremsstrahlung 

Consider bremsstrahlung of gluons with energy xE ^> from a high-energy particle 
of energy E and species s. The leading-log result for bremsstrahlung from this particle is 
(ignoring final state factors for the high-energy particle and gluon) 4 



dT s ^ gs _ a/i 2 ± P s ^ g (x) 



dx Any/2 x(l -x)E ' 



(2.12) 



where x is the bremsstrahlung gluon momentum fraction, P s ^ gs (x) is the usual vacuum 
splitting function, and, at leading-log order, 

til ~ [8x(l - x)E [\C K + (C s - \C A )x 2 + \C A {1 - x) 2 ] |(Q ±0 )} . (2.13) 



4 The NLL analysis will be based on Arnold and Dogan [11]. The leading-log result in (2.12) and (2.13) of the 
current paper roughly corresponds to eqs. (1.1), (3.1a-b), and (4.16) of Rcf. [11], with the correspondence 
that n± = todA^- The only difference is that we have written the answer in terms of the general q(Qj_o), 
whereas Ref. [11] specialized to the q± <C T limit, resulting in their Eq. (4.12). See Appendix D of the 
current paper for more detail, including an explanation of differences in notation. 
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Here, Q±o is any rough guess (1.5) of the total momentum transfer during the formation 
time, and ambiguities in that guess of 0(1) factors only affect the answer beyond leading- 
log order. The Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) splitting functions in 
(2.12) are 

*W*) = 11 + { \~ X?] , (2-14) 

P (X)-C. [1+3?+ (I-*) 4 ! (2l5) 

WO - C A ^— - } . (2.15) 

If we now follow the same steps as Arnold and Dogan [11], the NLL result corresponds 
to replacing (2.13) by 

III ~ [8x(l -x)E] 1/2 

1/2 



X 



where 



§C7 A $(^ ± ) + (C7- _I^ , (2.16) 



e = exp(2- 7E + f) (2.17) 

is the same constant as in Arnold and Dogan, and where (2.16) is an implicit equation 5 for 
/i±. Details are given in Appendix D. 

The assumption that goes into this result is that q(A) is proportional to In A at the scale 
of its arguments — that is, that dT e \/d 2 q± cc l/qj_ for q± of order the argument A of q. This 
assumption works for arguments in the range < A < T as well as the case T < A 
considered in this paper. As a result, if one uses the small q± form (1.4) of q(A), then the 
formula (2.16) reproduces the result found in Arnold and Dogan. For high energy particles 
with parametrically E ^> T/g 4 ln(l/g), however, one should instead use (2.6) for q. 

As discussed in Refs. [9, 11, 13], if the running of the coupling is included, the explicit 
factor of a s in (2.12) should plausibly be a s (fj,±). 

Our result (2.6) formally treated m D <C T. If one evaluated q(A) in a way that relaxed 
this assumption, Eqs. (2.12) and (2.16) could still be used to calculate the bremsstrahlung 
rate for high-energy particles. 

D. Result for Pair Production 

Pair production is the same as gluon bremsstrahlung except for (i) a change in group 
factors, and (ii) use of the corresponding vacuum splitting function 

P g ^(x) = N { t F [x 2 + (1 - xf}. (2.18) 



5 As discussed in Sec. Ill of Arnold and Dogan [11], one can choose to solve this equation iteratively, starting 
from some initial guess Qj_o and explicitly generating the first two terms of the expansion in inverse logs. 
However, as discussed in Ref. [11], the implicit form has the advantage that one does not have to generate 
an initial guess Q^o- 
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See, for example, the more symmetric presentation given in Arnold and Dogan [11]. The 
result (summed over quark flavors) is 



dY 



g->qq 



dx 



afi 2 ± P g ^ q (x) 
4ttV2x{1-x)E 



(2.19) 



with 



n\ ~ [8x(l - x)E] 



1/2 



X 



(C F - |C A ) $(^V) + |C A x 2 |(^) + §C A (1 - *) 2 J' 7 ' (2-20) 



III. dT el /d 2 q ± AND |(A) 
A. Strategy 

Our first goal is to determine how the differential elastic scattering rate dT e i/d 2 q± inter- 
polates at q± ~ T between the two limits shown in (1.3). For q± ^> tod, screening effects 
can be ignored, and the differential cross-section will have the form 

dT cl ^ _Cn_ x g^FjqJT) 



d 2 q± (2tt) 2 " q 4 ± 

for some function F(q±/T) with 

g 4 T 3 F(0) = g 2 Tm\ and ^ 4 T 3 F(oo) = (3.2) 

In the limit of weak coupling, once we know F(q±/T) we can then construct a formula valid 
to leading-order at all scales for q± as 

dT el _ C R x g*T*F{qJT) (3 3) 



d 2 q± (2tt) 2 ql(ql + m 2 D ) 

since did <C T. 

It will be convenient to write g 4 T 3 F in the form 



9 4 T 3 



g 4 T 3 F(q ± /T) = [~ h I + {qjT) + ~ i I_{qjT)Y—. (3.4) 

The functions I±(q±/T), to be determined, extrapolate between 

/±(0)=C±(2) and /±(oo) = C±(3). (3.5) 
The equivalence to (3.2) can be seen from Eqs. (2.3) for m D and Af. 
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B. Starting Point for dT e \/d 2 q± 



In general, the rate for a high-energy particle of energy E to scatter from the plasma is 
given by 

dY eX s f dq z y-^ _ f d 3 p 2 da c \ 



d 2 q A 



J -0)sJ2 d ^ J ^^ fs ^ P2) [1±/- 2 (P2-<?)L (3-6) 

S'2 



where z is the direction of motion of the high-energy particles, p 2 is the momentum of a 
particle in the plasma of species s 2 , a e \ is the cross-section for scattering from that plasma 
particle, /(P2) is a Bose or Fermi distribution that accounts for the probability of encoun- 
tering the plasma particle, and 1 ± f(p2—q) is a final-state Bose enhancement or Fermi 
blocking factor for the plasma particle after transferring momentum q to the high-energy 
particle. We assume that E ^> T and so do not need to include any final state factor for 
the high energy particle. In (3.6), v S2 is the number of spin and flavor degrees of freedom 
for species s 2 (2 for gluons, 4iVf for the sum of quarks and anti-quarks). 

For high-energy particles (in this case, E ^> mo), elastic scattering from the plasma is 
dominated by t-channel gluon exchange. The infrared behavior of t-channel gluon exchange 
is cut off in the infrared at qj_ ~ by the effects of Debye screening and related phenomena 
in the plasma. To leading order in the weak-coupling limit, the problem is simplified by the 
fact that we can ignore screening effects when investigating q± ~ T since T 3> m D . To 
leading order, the differential elastic scattering rate for q± 3> m D is then 



dT e i _ f dq z duj ___ f d 3 p 2 C R t R2 g 4 



A' 2 



d 2 q± J (2vr) 4 ^ 62 J (2tt) 3 (2p) 2 2p 2 2\p 2 - q 



Q 1 



f S2 (p 2 ) [1 ± f S2 (p 2 - q)} 



X27r5(uj-q z )27i5(uj+\p 2 -q\-p 2 ) (3.7) 



for a massless high-energy particle with momentum p in the z direction. We use capital let- 
ters for 4-momenta, with P = (p,p) and Q = (cu, q). The factors of (2p)~ 1 , (2p) _1 , (2p 2 )~ 1 , 
and {2\p 2 — q\)~ l are the usual initial and final state relativistic phase space normalizations, 
where we've taken the high energy limit E 3> q. The 5(u + \p 2 — q\ — p 2 ) is energy con- 
servation for the (massless) plasma particle. The other 5-function is energy conservation 
$(\P + Q\ ~ 00 ~ P) f° r the incident high-energy particle, again taking the high-energy limit 
p ^> q. The two gluon vertices in the amplitude are 2gP^ and 2g P!^ (times color generators), 
regardless of the types of particles colliding. The appearance of this universal form can be 
understood as a consequence that, in the high-energy limit E ^> T, the upper limit (1.5) 
for the range of individual q± transfers dominating bremsstrahlung is small compared to the 
center-of-mass energy ~ (ET) 1 ^ 2 for a collision with a plasma particle. In the center-of-mass 
frame, the t-channel gluon is soft compared to either particle involved in the elastic collision, 
and so we may use the universal form that gluon-particle vertices take in the soft gluon limit. 
Note that P^P 2 in (3.7) could have equally well been written using the final plasma particle 
momentum P 2 — Q as P^{P 2 — QY, due to the 5(cu — q z ). 
Performing the 00 integration in (3.7), 

o?r e i C R ^ 4 f dq z f d 3 p 2 (p 2 - p 2 . 



d 2 q ± (2tt)M 



E- . 4 f dq z f d 3 p 2 (p 2 - p 2z ) r , 

^ j ^ J J2nf ^p^ UM ^ ± ^ ~ ^ 

x27r5(q z + \p 2 -q\-p 2 ). (3.8) 
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In thermal equilibrium at zero chemical potential, the distributions f S2 are all the same for 
massless bosons, and also all the same for massless fermions, and we can rewrite dT e i/d 2 q± 
as 

with 

x2n5(q z +\p 2 -q\-p 2 ). (3.10) 

C. Recasting dT e \/d 2 q± as a double sum 

Now rewrite the equilibrium Bose and Fermi distribution functions as sums of exponen- 
tials, 6 



f±(p) = -p—r = ^(ir-V^, (3.11) 

m=l 

oo 

1 ± f ± (p) = J2(±) n e~ nPp , (3-12) 



ra=0 



with P=l/T. Then 



with 



J ± (g ± /T) = J2 E( ± ) m+n_1/ -(^/ T ) ( 3 - 13 ) 



m=l ri=0 



T 3 y 2tt (27r) 3 p 2 |P2 - <?| 
This integral is evaluated in Appendix A and yields 

where 1^(2) is the modified Bessel function of the second kind. The case n = gives the 
n — > limit of the above formula, 

4o = ^ ■ (3.16) 

m 3 

In the gj_ — > limit, J mn — > l/(m + n) 3 , and so 

/ ± (0) = VV^ ^= V P ^ = C ± (2). (3.17) 

±v 1 A^A^( m + n )3 (m + n) 2 v ' v ' 

m=l n=0 m+n=l 



6 The analysis of Rcf. [11] could also be applied to non-equilibrium isotropic situations. Here, however, we 
are specializing to equilibrium distributions. 
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In the opposite limit of q± — > oo, only the n = terms survive in the double sum (3.13), 
giving 

oo 

I±(oc) = ^(i)"- 1 /^ = C±(3). (3.18) 

m=l 

These two limits are in accord with (3.5). We will later find it useful to extract the n = 
contribution from the general case, writing 

I±(q±/T) = J±(oo) + AI±(qjT) = C±(3) + AI±(q ± /T) (3.19a) 

with 

oo 

AI±(qjT) = (±)™I mn (qjT). (3.19b) 

rrt,n=l 

If one wished to evaluate the functions I±(q±/T) numerically, the expansion (3.19b) 
converges rapidly for q± > T. For the case of small q±, we show in Appendix B that the 
expansion is 7 

I + (Q) = C+(2) [l-Jfl + (i + ^)Q 2 + 0(Q 3 )] , (3.20) 
I-(Q) = C-(2) [1 + (£ - i^)2 2 + ^(Q 3 )] • (3.21) 



D. Integrating to get q 

Using (3.3) and (3.4), our leading-order formula for the differential elastic scattering rate 



is 

dT el ^ C R g^[~ h I + {q ± /T) + ~ i I-{q±/T)\ 



d?q ± (2vr) 2 n 2 ql(ql + rnl) 

So, to compute the integral (1.1) that gives <f(A), we turn to evaluating the integrals 



(3.22) 



X± (A)./ ^^X. (3.23) 
J q± <A (2tt) 2 (ql + m 2 D ) 

In the weak coupling limit, one can choose a momentum scale A between ttid and T, with 
m D C A < T. Then one may split the integral at A into separate pieces for which q± <^ T 
or q± ^> mo approximations may be made: 

J±(A) ^/ g$. U°) f p^'^JH.. (3.24) 

In Appendix C, we then evaluate these integrals using the double sum (3.19) for the last 
one. To leading order, we obtain the result (2.6) with 

2 ^ \ i (m + nY ^ K ' (m + n) 3 v ' 

m,n=l v ' m,n=l 



7 A numerically modest formula which reproduces I+{Q) to a few tenths of a percent is to (i) for Q < 3.2 
use (3.20) plus —0.0062 Q 3 inside the square brackets, and (ii) for Q > 3.2 use the m + n < 3 terms of 
(3.19). For I-(Q) do the same, but with -0.0031 Q 3 added inside the square brackets of (3.21). 
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This can be written as <7+ = — <9 a T(a, 0, 3)| where 

x 1 

T(a,b,c)= V j- (3.26) 

^ m a n b (m + n) c 

m,n=l v ' 

is the Tornheim zeta function [16], and similarly <r_ can be defined in terms of its general- 
ization. However, this does not seem to give any information that is more useful than the 
sum itself. By letting k = m + n and then summing over n for fixed k, one obtains the 
single-sum formula (2.6d) for a±. 
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APPENDIX A: I mn 

In the integral (3.14) defining I mn , use the 5-function to rewrite 

(P2 - P2z) 2 = (P2 - P2z) [\P2 ~q\- (P2z ~ Qz)] • (Al) 

Then rewrite the 5-function as 

/+oo 
dAe afe+| P2 -q|- P2 )_ ( A2 ) 
-oo 

The p 2 integral now takes the form of a convolution of functions W m and W* defined by 

Wm (p 2 , A) = (P2 ~ P2z) e -( m/3+iA ^. (A3) 

P2 

To turn the convolution into a simple product, we Fourier transform from p 2 to r: 

Imn(qJT) = ^J^±J d \ J d\ W m (r, A) W* n {r, A) e~^\ (A4) 
The q z integral then gives 8(X — z), which can be used to do the A integral: 

2 r 

I mn {qjT) ee ^ y dh W m (r, z) W:(r, z) e^ r \ (AS) 

The Fourier transform of (A3) evaluated at A = z is 

T 3 m 

W m {r,z) = — — — (A6) 

[m 1 + 2imzT + (r^T)^J^ 

Now do the z integration in (A5) by closing the contour in the upper-half complex plane 
and picking up the double pole there: 

U«i/^ = —! d 2 r± , - , 3[ ( rt^ 13 e ~ iq± ' r± - ( A7 ) 
7r y (m + n) i [(r ± T) 2 + mn]" 3 

Finally, performing the r± integral gives (3.15). 
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APPENDIX B: SMALL q ± EXPANSION OF I± 



It is possible to find the small Q expansion of I±(Q) starting directly from the integral 
formula (3.10). For the bosonic case, at least, we find it easier to instead start from the 
double sum formula derived in Sec. Ill C. 



1. The 0(q±) piece of 1+ 

Start from the double sum of (3.13) and (3.15) and subtract off the q± = piece: 

m,n=l ^ ' 

For small Q = q±/T, this sum is dominated by large m and n, and so we can replaced the 
sum by an integral: 



3C 



51+ ~ / rfm rfn 7 V \ q ; i . B2 

7o m + n 3 



Change integration variable from n to x = Q^mn, and then do the m integral to get 

2 



5I + ~^Q I dx 
o 



*2(*)- ; , 



-ge. (B3) 



This gives the O(Q) term in (3.20). 
2. The 0(g5_) piece of I + 

The difference between the double sum and the integral approximation made above is 

ei + =(±-r *» dn ) agwgvp - a . (B4) 

io / (m + n 3 

Now rewrite 

££-// = (£-/)/+/(£-/) + (£-/)(£-/> (b 5 ) 

and correspondingly (using the m <-> n symmetry of the sum for J + ) 

5 2 / + = 5 2a / + + 5 2a / + + 5 2b /+. (B6) 
To evaluate S 2a I + , again change integration variable from n to x = Q^/mn to get 

I oo 

5 2a / + = Q 4 ( V'_ / dm)F(m;Q) (B7) 



£-/ *») 

m=l ^ / 
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where 

, f°° dx x[\x 2 K 2 (x) - 1] 
For small Q, F(m; Q) is a slowly varying function of m, which in general means that 

V-/ dm )F(m; Q) ~ -±F(0; Q). (B9) 

m=l ^ ' 

In our case, F(0; Q) should be understood to be the m — * limit — 1/8 Q 2 of (B8), which 
then gives 

5 2 *I + = ±Q 2 . (BIO) 

The other term 5 2h I + in (B6) will be turn out to be dominated by m and n with ran <C 
I/O 2 in the small Q limit. Making this small Q approximation to the argument of K 2 , we 
get 

( °° poo \ / °° /"OO \ 

V-/ ^m) V-/ dn)G(m,n) (Bll) 
m=l ^° / \n=l ^° ' 

with 

g(m,n)= - (B12) 
4(m + n) 3 

Now rewrite 



°° /"OO \ / /"OO \ 

— / dm ] I — / J G(m, n) = 

m=l ^ / V n=1 -A) / 

\ G(m, n)- dm 1 G(m', n) - dri G(m, ri) + dm' dri G(m', ri) 

m n= i I J m— 1 J n— 1 Jm—1 J n— 1 

(B13) 

Explicit evaluation of the integrals using (B12) yields 



c-9h ~r ^9 \ (m — 3m n — 3mn + n ) + Amn 

S Zb L = Q 2 > -i- — — — - , (B14) 

„ x 8(m + n) 3 (m + n- l) 2 (m + n-2) v ; 

where the summand for m = n = 1 should be treated as the limiting value — ^. Defining 
= m + n, (B14) can be rewritten 

+ fe^ 8*3(ifc- 2) • (B15) 

Doing the n sum first (and treating the k — 2 case separately, which does not fit the pattern 
ofA;>2): 



5 2b I + = -Q ^ 

k=3 



^ + ^2ikHJ-T)\ = - Q '^-^ (B16 » 



Putting (BIO) and (B16) into (B6) then yields the 0{Q 2 ) term in (3.20). 
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3. The 0{q 2 ± ) piece of i_ 

The case of J_ is much simpler because there is no 0(Q) term in the expansion. A quick 
but non-rigorous way to obtain the answer is to naively expand the summand of (3.13) and 
(3.15) in powers of Q: 



^- = -Q 2 E (-) 



m+n-l mn 



Aim + n) 3 

m,n=l v ' 

oo k— 1 /, \ 

fc=2 n=l 

= -^2 2 E(-) fe - 1 (l-^)- (BIT) 

k=2 V ' 

If one then interprets XX ~ ) fe_1 as being C-(0) = §, then 

S 2 I- = -UU0)-C-mQ 2 , (B18) 

where C_(2) = 7r 2 /12. This gives the result (3.21) quoted previously. 

A more reliable way to make the same calculation is to start with the convergent, un- 
expanded sum analogous to (Bl), 

oo 

j_(C) - Mo) = £ (-) m+n-1 #K™; 2), (Bi9) 

m,ri=l 

\lQ 2 mnK 2 (Qx/rrm) - ll 
H(m, n; Q) = ^ 2 ^ 3 ; ± . B20 

Now block the sum into 2x2 blocks as 
J_(Q)-J_(0) = £ [-#(m,n;Q) + #(m+l,n;Q) 

m,n odd 

+ if(m,ra+l; Q) - H(m+l,n+l; Q)]. (B21) 

At this point, one can safely expand the summand to order Q 2 . Then change summation 
variables from m to k = m + n, and then sum over first n and then k. The final result is 
(B18). 



APPENDIX C: X± 

In this appendix, we evaluate the integrals in (3.24). The non-trivial integral is the second 



one: 



T(2 ) _ / dV I±(qJT) Moo) /A\ f A/ ± (g) 

" WA (27T) 2 " £ 27T U; M< S</ 3A (2VT) 2 Q 2 1 J 
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where the last equality uses (3.19) and switches to the dimensionless integration variable 
Q = q±/T. Because Ai± falls off for II — > oo, we can drop the UV regularization Q < (3A 
in the last integral. 

It's useful to now change infrared regularization by inserting an initially unnecessary 
factor of Q 2t with the limit e — > + taken at the end of the day. Then rewrite the above as 



T ( 2 ) _ / ± (oo) (A\ f d?Q AZ±(g) f d?Q AZ±(g) 

J± -^r ln UJ + y w Q^"y s <,A(^F (C2) 

In the last term, we can replace Ai±(g) by A/±(0), giving 

(C3) 



/ 

Jo, 



d 2 Q AI ± (Q) _ _ (/3A) 2 < _ A/ ± (0) 



l + 21n(A) + 0(r) 



' 2 </3A (2vr) 2 Q 2(1 ~ e) ±V 7 47re " 4tt 

For the other integral in (C2), we use the double sum formula of (3.19b) and (3.15) for Ai±: 
f d*Q AZ±(g) _ ^ t mn f d 2 Q , x 

OO / \ g /*oo 

; 4tt (m + ra) 3 7 

m,n=l 



The last integral gives 2 2e T(e) T(2 + e). Expanding the result in e, and noting that (3.17) 
and (3.18) give 

00 /_i_\m+n— 1 

E 7 - v = C±(2) - C ± (3) = A/ ± (0), (C5) 
' (m + nr 

m,n=l 

we get 

d 2 Q AZ±(g) AZ±(0) /l \ <t ± 

with a± defined as in (3.25). 

Evaluating J ± = j£ } + j£ 2) of (3.24) by combining (C2), (C3), (C6) and 



± " X ± <a (2vr) 2 (g 2 + m 2 ) 2vr UJ 1 J 

then produces the final result (2.6). 

APPENDIX D: THE NLL CALCULATION 

For simplicity of presentation, we will focus just on bremsstrahlung in this appendix. 
Arnold and Dogan [11] computed the NLL bremsstrahlung rate using the q± <C T limit in 
(1.3) for dT e \/d 2 q±. In their notation, they referred to A(q±) instead of dT e \/d 2 q±, and the 
translation 8 is 

wrU AM < (D1) 



See Appendix A of Ref. [13]. 
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9 2 I |% A(q ± ) ql (D2) 



(2vr) 2 

The derivation of Arnold and Dogan is fairly easy to generalize if we just make these replace- 
ments. If one is confident enough, one can just write our generalization (2.16) of Arnold and 
Dogan by inspection by recasting their result in terms of the A < T version of q and then 
assuming the formula works for A>T. This works because Arnold and Dogan's constant £ 
was generated by the large-g_i_ contributions to the calculation, where A(q±) is proportional 
to \jq\_ in either case. It is the large q± part of the calculation that is affected by the 
NLL calculation: the qj_ <C Q± pieces just come from the leading-order calculation, which 
is proportional to q(Q±o). 

Readers may find the above argument obscure, or an explicit calculation reassuring, and 
so we will also indicate how to get the same result by modifying the calculations of Ref. 
[11]. We will not reproduce the entire derivation of Ref. [11] but will just indicate which 
equations are modified. 

Eq. (4.14) of Ref. [11] for H 2 , which determines the leading log result, becomes 

H 2 = {2g 2 \p'kp\ [\C A p 12 + (C s - \C A )k 2 + \C A p 2 ] q-(Q ±0 )} 1/2 , (D3) 

where p' — E is the initial high-energy parton energy and k = xE and p — (1 — x)E are 
the energies of the two partons it splits into. Combining this formula for H with Eqs. (1.1), 
(4.1), and (4.15) of Ref. [11] gives Eqs. (2.12) and (2.13) of the current paper. Roughly 
speaking, H gives the scale of the total momentum transfer to a high-energy particle during 
the formation time as 

Qx ~ ^, (D4) 

where Xi — 1, x, or 1 — x depending on which particle one is focusing on. 

In evaluating the NLL correction, the only remaining step that is different is the evaluation 
of the integral 

Hk2) = g2 1 (0 A{q±) Fo(h) ' [Fo(h) - Fo(h + Kq±)] (D5) 

defined in Eq. (4.27) of Ref. [11], with 

F (h) = iAp'kp 



exp ( -e W4 ^) 



and k = p', p, or k. Here also, Arnold and Dogan used the small-g_i_ formula for A(q±_), and 
we now want to generalize. Introduce a cut-off A <^ Q± such that the differential elastic 
scattering rate behaves like 

§^ = ^§^^^70^ ^rA< q ± <Q ± (D7) 
d 2 q± d^qt [2n) 2 q\ 

for some constant c, up to corrections higher-order in g. We won't need the specific value, 
but, referring to (1.3), c = g A M for the case Q± ^> T of interest in this paper, or c = g 2 Tm^ 
for the case <^ Q± <^ T analyzed previously in Ref. [11]. Now rewrite 

I 2 (k 2 )^I 2< (k 2 ) + I 2> (k 2 ), (D8) 
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with 



'2<' 



/ 2 >(«: 2 ) = 



» 5± <A A± 

(i 2 /i /" oPg^ c 



(2tt) 2 7 A<?± (2vr) 2 gl 



- T F (/i)-[Fo(/ l )-F (/i + ^ ± )]. 



(D9) 
(D10) 



In the integration (D9) for 7 2 <, we have Kq± <C kQ± ~ and so one can expand the 
difference of F 's in a Taylor expansion: 



! 2< 



> 2 ) 



r d 2 h r 



d 2 q± 



dT el {q. 



± 



9±<A 
2/ 



giFo(h) • V 2 F (h) 



--k 2 I(A) (p» 2 ^- 

7T _H Z 



Now turn to I 2> . It is convenient to rewrite 

h> = 

with 

c 



1 2> ' 1 2> 



(Dll) 



(D12) 



) 2 



' 2 > J (27r) 2 7 (27T 

/• d 2 /i r d 2 q ± 
J (2?i 



gl % ^ A) "^i + M 2 ) 



1 2> 



(2tt) 2 g 2 (g 2 + M 2 ) 



F (h) • [F (h) - F (h + /cgj] , 

(DI3) 

F (h) ■ [F (h) - F (h + Kq ± )} . (D14) 



Here 9(z) is the step function and M < A is an arbitrary scale. In the first integral, we can 
again treat Kq± as small and Taylor expand the difference in F 's, giving 



J 2> — 



d 2 qj_ ( c 



0(q± - A) - 



ij (27T) 2 V?l 

T x 27 111 (t) x 



ql(ql + M 2 ) 

in/A 



(2tt)= 



Fo(fc) • V 2 F (h) 



(DI5) 



is proportional to the hi^ 2 ) integral evaluated by Arnold and Dogan, with m 2 ^ replace 
by M 2 in the denominator, and the overall normalization g 2 Tm^ replaced by c. So we can 
take over the result from Eqs. (4.31) and (4.32) of Ref. [11], 



J 2> — 



(p'kp) 2 c 

7l 2 M 2 



(2 - 7e - ln« K )w K 



with 



2,„2 



u K = e"/ 4 



Combining the various pieces above, 
2K 2 e in ^(p'kp) 2 



h(K 2 ) ~ 



vriJ 2 



<< A > " s 



2# 2 



-2 + 7E + ln(e-/ 4 ^) 



(D16) 
(D17) 

(D18) 
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which replaces Eq. (4.32) of Ref. [11]. Now take the real part of I 2 and note that (1.1) and 
(D7) implies that 

|(A') ~ |(A) + ± In (£) (D19) 

for A' and A both in the region covered by (D7). Using the definition (2.17) of £, we then 
get 

which replaces Eq. (4.33) of Ref. [11]. Eq. (4.35) of that reference then becomes 
Re(S, F s ) = Re(S, F ) 

\C,P- q + (C s - \C K )1? q ( ^f) + q (^) 



Re/ 2 (« 2 ) ~ ~ V :U g \/^H, (D20) 



X 2 i 1 + ±C A p' 2 q(Q ±0 ) + (C s - \C A )V |(Q_l ) + |C-AP 2 |(Q±o) 



(D21) 



where the notation (S,F S ) is defined in Ref. [11]. Following the same steps as Arnold and 
Dogan then produces our result (2.16), where our fi± corresponds to their mj)fi±. 
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